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ABSTRACT 

We study the effects of strong lensing on tlie observed number counts of mm sources using a ray 
tracing simulation and two number count models of unlensed sources. We employ a quantitative 
treatment of maximum attainable magnification factor depending on the physical size of the sources, 
also accounting for effects of lens halo ellipticity. We calculate predicted number counts and redshift 
distributions of mm galaxies including the effects of strong lensing and compare with the recent source 
count measurements of the South Pole Telescope (SPT). The predictions have large uncertaities, 
especially the details of the mass distribution in lens galaxies and the finite extent of sources, but the 
SPT observations are in good agreement with predictions. The sources detected by SPT are predicted 
to largely consist of strongly lensed galaxies at z > 2. The typical magnifications of these sources 
strongly depends on both the assumed unlensed source counts and the fiux of the observed sources. 
Subject headings: galaxies: high-redshift, Galaxy: formation, Galaxy: structure, gravitational lensing: 
strong, methods: statistical, submillimeter: galaxies 



1. INTRODUCTION 

Many galaxies at redshifts z^2 — 5 have been found to 
be undergoing large amounts of star formation, leading 
to a population of distant galaxies with large amounts 
of warm dust that can be observed at in frared, submm, 
and mm wavelengths (jBlain et al.l 120021 ). Star- forming 
rates are found to be o ften in excess of 1000 Mq/jt 
(jMichalowski et aLll2010[ ). contributing a significant frac- 
tion of the total cosmic star formati on at these redshifts . 

Surveys at submm wavelengths (jCoppin et al.l l2006( ) 
covering smaller areas at high sensitivity have estab- 
lished the existence of a population of dust y star-forming 
galaxies at high redshift (|Blain et al.ll2002l ). showing that 
the number of sources as a function of flux (the lu- 
minosity function) is steeply declining at high fiuxes. 
As a result, gravitational lensing is expected t o sig- 
nifica ntly modify the observed number counts (|Blainl 
|1998[) . an eff ect known as "mag nification" or "amplifi- 
cation" bias ()Turner et al.lll984) . Recent mm- wave sur- 
veys like the South Pole Telescope (SPT; Vieira et al 
2010) are now covering enough area to accumulate sta- 
tistically significant numbers of highly luminous distant 
galaxies, providing an opportunity to compile large sam- 
ples o f strong gravitational lenses. Recent theoretical 
work (iNegrello et al.|[20q7t iFedeli fc Berciano Albal [200l 
iLima et al.ll2009l . 120101 iJain fc Limall2010( ) hasdemon- 
strated that gravitational lensing is likely an important 
contributor to the galaxy counts observed by large scale 
mm- wave surveys. In addition , evidence is now emerg- 
ing from Herschel observations (jFraver et al.ll2010() that a 
large fraction of the brightest high-redshift dusty galaxies 
are indeed strongly lensed. 

Much remains unknown about massively star-forming 
galaxies at high redshift. The redshift distribution as a 
function of flux is only roughly understood, and differ- 
ent models have very different input physics. For exam- 
ple, the Durham semi-analytic model of galaxy forma- 
tion requires a top- heavy IMF to explain this population 
(jLacev et al.ll2010[ ). Strong lensing of these sources al- 



lows a magnified view, making multi-wavelength follow- 
up easier, as the sources are brighter. 

In this work we calculate the expected number of 
strongly lensed galaxies in flux-limited mm- wave surveys, 
paying particular attention to the expected redshift dis- 
tribution of the sources and lenses and the effect of finite 
source effects. 

2. OVERVIEW OF CALCULATION OF NUMBER OF 
LENSED SOURCES 

Determining the expected number of galaxies discov- 
ered in mm- wave surveys is complicated by at least four 
major uncertainties: 

• the statistics of the source population (uncertain 
number cotmts, uncertain redshift distribution) 

• the properties of the source population (uncer- 
tain spectral energy distributions, uncertain angu- 
lar sizes) 

• the statistics of the lens population (number counts 
as a function of mass and redshift) 

• the properties of the lens galaxies (internal mass 
profiles and ellipticities). 

To investigate uncertainties in the source population, 
we use two independent unlensed source count predic- 
tions for SPT measurements at 220 GHz (1.4 mm). In 
particular, the redshift distribution is expected to play 
a key role in determining lensing efficiencies, so the two 
different models are intended to provide an estimate of 
the redshift importance. Our first model, henceforth 
called t he Durham model, based on the models devel- 
oped in iBaugh et al.l (|2005[ ) (Lacey et al. private com- 
munication), is the result of semi-analytic modeling of 
galaxy formation. The second model consid ered, referred 
to as the UBC model ()Marsden et al.|[20Tol) . is obtained 
through backward evolution models of the local Universe. 

To model the lens population we follow 'Perrotta et alj 
|2002), using a Press-Schechter (,Prcss fc Schechtcr„197¥ 
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approach to the lens halo distribution as a function of 
mass and redshift. We use the Sheth & Tormen (1999) 
redshift-dependent mass function to model the number of 
halos of a given mass and redshift for our lens population. 

For the internal mass distribution we assume that the 
region where the majority of strong lensing occurs can be 
modeled as an elliptical mass profile with a 3D density 
profile that falls as This is an exce llent approxima- 

tion for galaxies (jKoopmans et al]|2009D , while it is likely 
not a sufficiently complex model to capture the lens- 
ing p roperties of massive galaxy clusters (Richar d et alJ 
1201(1 . 

We use ray-tracing simulations to explore the impact of 
lens ellipticity and finite source sizes, assuming constant 
values of ellipticity and source sizes and exploring the 
impact of different assumed values. 

In all that follows, we assume as our fiducial cos- 
mology a spatially flat universe with rim = 0.222, 
Ho = 71.0kms-iMpc-i, Qb = 0.0449, = 0.963 and 
as ^ 0.801. 

3. UNLENSED MM-WAVE NUMBER COUNT PREDICTIONS 

The assumed unlensed source count models (the 
Durham and UBC models) have not been calibrated at 
mm wavelengths; small differences in parameters such as 
dust emissivity that are not large effects at submm wave- 
lengths could lead to large misestimates at mm wave- 
lengths. 

A simple check is to verify that the models produce a 
reasonable amount of noise power in mm- wave maps ; this 
has been measured in SPT data in lHall et al.l (|2009| ). We 
computed the angular noise power, defined for randomly 
distributed point sources as 

ci=j^ S'^d{lnS) (1) 

with a flux cutoff of 17 mJy (|Han et al.l 120091 ) for SPT 
at 220 GHz (1.4 mm). The majority of the noise power 
comes from sources well below the SPT sensitivity limit 
for detecting individual sources. 

The noise power from the UBC model is in excel- 
lent agreement with the measurements, and did not re- 
quire any corrections. The noise power derived from the 
Durham model, around 60 Jy^/sr, i s more than a factor 
of 3 too high compared to lHall et al.l (|2009[ ). We scale the 
flux of each object by a constant amount to match the 
measured noise power. In Figure[l]it can be seen that the 
Durham and UBC models are forced by this constraint 
on the total power to have comparable number counts at 
fluxes of a few mJy. 

After applying a constant flux correction to match the 
observed noise power, the Durham model still had a 
problem with the properties of low-redshift dusty galax- 
ies, in that the predicted number of low-redshift galax- 
ies in the SPT sample was too high. In particular, the 
brightest of the low-z population in the Durham model 
(observed by IRAS) are predicted to be bright enough 
in the SPT maps to be found as sources. This is not 
the case, as evidenced by the small fraction of SPT- 
discover ed galaxies that w ere also observed to be IRAS 
sources (jVieira et al.ll2010l P1. As gravitational lensing is 

^ http:/ /polo. uchicago.edu/public/data/vieira09/indox. html 



more efficient for high redshift sources, we elected to sim- 
ply further suppress the flux of low-z galaxies {z < 0.2) 
to make the low-z galaxy counts agree with the num- 
ber of SPT sources found to coincide with IRAS sources. 
While not rigorously justifiable, the main point of using 
the Durham model was to get a plausible redshift distri- 
bution of sources at high- redshift. In order to avoid mis- 
interpretations caused by the low-z population we only 
apply our lensing model to Durham counts with z > 0.2 
eliminating the low redshift counts and compare the re- 
sults to SPT's IRAS removed counts. 
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Fig. 1. — Unlensed number count predictions for SPT sources 
at 220 GHz. The thin dash-dotted blue line shows the original 
Durham counts while the thin dashed blue curve shows the same 
counts with ^ > 0.2. The thin solid blue curve is t he Durham mode l 
corrected to match the measured noise power HHall et al.ll2009l V 
The thick black dash-dotted line shows the total UBC counts and 
the thick solid black line is UBC with z > 0.2 for comparison pur- 
poses. The data points are 220 GHz (1.4 mm) SP T dusty galaxies 
after removing known IRAS sources jVieira et al.ii2010.) . 



4. GRAVITATIONAL LENSING THEORY 

For the details of gravitatio nal lensing theory we re- 
fer th e reader to a review by IBartelmann fc Schneiden 
(|2001f) . Here we only briefiy state a few lensing quanti- 
ties that are used in this work. The lens equation that 
we solve numerically is written as 

/3 = 9- ^a{Dj) = 9- a{9) (2) 

where Dj,^ Dg, and D^g are the angular diameter dis- 
tances of deflector (d) and source (s) and /3 is the ob- 
served position of a point source at 9 deflected by an 
angle a. The total deflection from an ensemble of point 
masses at a single lens plane is given by 

d{e) ^ - f d?9'K{9') (3) 
TT J \0-9'\'^ 
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The dimensionless surface mass density k is 
'^(^) — — ^ where Scr = * 



4nG DdDds 



(4) 



is the critical surface mass density and is the 2D 

projected mass density of the lens. 



5. LENS MODELING AND CALCULATION DETAILS 
5.1. Ray-tracing 
The lens equation (Equation [5]) is an implicit equa- 
tion, in that the image position 6 is needed to evaluate 
the deflection d{6). For a given source position (3 it is 
difficult to solve this implicit equation to obtain 9. On 
the other hand, if 6 is known /3 can be easily evaluated. 
This constitutes the basis of ray-tracing simulations of 
gravitational lensing: start with an array of image posi- 
tions and determine the source positions from which they 
emerged. 

Simple halo profiles, like the singular isothermal sphere 
(SIS), can often produce analytical cross-sections. More 
complex mass profiles often have no simple analytical so- 
lution and the cross-section has to be numerically com- 
puted using a ray tracing simulation. In addition to the 
ability to compute lensing quantities for any arbitrary 
mass configuration, ray-tracing simulations have the ad- 
vantage that they can easily include finite source effects 
which are not properly modeled in analytical solutions. 

Our simulation makes surface density maps as a matrix 
of 400 X 400 elements a nd solve s the l ens equation using 
the method described in lKeetonI ()2001[ ). The image plane 
is divided to 180 x 180 squares each of which is divided 
into two triangles and the corresponding positions of each 
vertex are found in the source plane. The magnification 
is computed by defining a grid in the source plane and 
calculating the image-plane area of triangles that have 
been mapped into each source position. 

5.2. Halo Mass Profile and Ellipticity 

We assume that the lens profiles in the region of inter- 
est are well-approximated by singular isothermal profiles, 
where the three dimensional density profile for a spheri- 
cal profile is 

- ^ (5) 

where cr^ is the line of sight velocity dispersion of the 
stars in the galactic disk or the galaxies in a galaxy clus- 
ter. It is well known that this is not a good approx- 
imation to halos pr oduced in dark matter simulations 
([Navarro et al.|[T997l ). but the region that dominates the 
strong lensing properties is typically dominated by bary- 
onic proces ses and is empir ically found to be close to 
isothermal (jKoopmans et alll2009[) . 

The dependence of the velocity disp ersion on the red- 
shift a nd mass of the halo is given by iBrvan fc Normaiil 
IpM) as 



a, = /.Ml/3 



g^(z)A(z)G'' 
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where 

A(z) = IStt^ + 82{Qrn - 1) - 39(an - 1)^ (7) 

and fa- is a scaling parameter used to match the nor- 
malization from simulations. While the density profile 
may indeed scale roughly as 1/r^ in the central regions 
of interest, the relationship between the normalization 
(i.e., the velocity dispersion) and total halo mass is not 
empirically calibrated and is a potential large source of 
systematic uncertainty in this analysis. 

Integrating p{r) along the line-of-sight produces the 
projected surface mass density 



2G^ 



(8) 



The corresponding dimensionless surface mass density is 



- To 



(9) 



where we have defined the Einstein deflection angle as 



= 47r 



/ o-iy \ ^ Dds 
\ c) 



(10) 



The magnification for point sources lensed by this mass 
profile is analytic, with ^{9) = \9\/{\9\ — 9e)- We calcu- 
late magnification for extended sources as the ratio of the 
area of combined images to the source area. This allows 
easy exploration of finite source effects and non-trivial 
lens distributions. 

The integral lensing cross-section a{fi > fimin) is the 
area on the source plane inside which the magnification of 
a source is equal or larger than ^mm- Throughout this 
work > fimin) and a{^) are used interchangeably. 
The cross-section for an SIS halo for ^min > 2 is given 
analytically by Perrotta et al. (2002) and at large mag- 
nifications our numerical results agree with the analytic 
form to better than 2.5%. 

The morphology o f the galaxies from observations 

vans fc Bridlell200"9f) and th e shape of dark matt er ha- 
los from N-body simulations (jLudlow et al.|[20Tol) both 
indicate that a considerable amount of ellipticity is 
prese nt in the lensing halos. Many studies previously 
fe.g.. iMeneghetti et al.l (120051) 1 have introduced elliptic- 
ity in the projected two-dimensional lensing potential ip. 
For analytical studies this has the advantage that the 
second derivatives of the potential directly give the lens- 
ing quantities and lead to simple analytical expressions. 
However for high values of ellipticity this implies dumb- 
bell shape density profiles which are unrealistic. Simu- 
lations and observations both indicate ellipticity in the 
distribution of mass, rather than in the potential result- 
ing from it. In our ray-tracing approach it is simple to 
introduce ellipticity in the actual surface mass density. 
At each point defined by x and y on the two-dimensional 
plane with radius r = x"^ + y'^ we compute the SIS or 
NEW density using radius re defined as 



+ y^(i-e«) 



(11) 



The axis ratios are now given by 6/a = 1 — e„. 

The impact of ellipticity can be clearly seen in Figure 
H EUipticity increases the area of the source plane where 



4 



Hezaveh, Holder 



strong lensing can occur. However, in the same figure we 
can see the effects of finite sources. The regions of high 
magnification become narrower as ellipticity increases, so 
larger sources will tend to have smaller peak magnifica- 
tions. 

The interplay between ellipticity and source size is 
shown in Figure [31 Ellipticity generally leads to higher 
magnifications, but finite source sizes become more im- 
portant for higher ellipticities. This discussion has as- 
sumed a fixed Einstein radius. The relevant factor is the 
ratio of the source size to the Einstein radius; a given 
source size that may be a problem for a galaxy-scale lens 
would behave like a point source for the purposes of lens- 
ing by a galaxy cluster. 

At intermediate magnifications, finite source effects 
tend to increase the lensing cross-section. This occurs 
because sources are now occupying a larger fraction of 
the source plane and are more likely to have a part of 
the emitting region be in the strong lensing regime. 
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Fig. 2. — Magnification maps of two halos with identical total 
mass (M = 10 Afo) and different ellipticities (e^ = 0.2 for tfie 
left panel and = 0.3 for tfie rigfit panel). Tfie extended circular 
source (drawn as a black circle) is identical in both plots making all 
magnification differences due to ellipticity. The left panel contains 
large magnifications in excess of /i = 100 which are damped in 
the right panel. But intermediate magnifications are extended to 
larger regions in the right panel. 



5.3. Probability function 

We use the halo mass function of iSheth &: TormenI 
(|1999f) to determine the number of halos of a given mass 
at each redshift, sum the source plane cross-sections of 
all lenses between the observer and source, and divide it 
by the total area of the sphere centered at the observer 
with a radius at the source redshift. This is the prob- 
ability that a source at that redshift is magnified by at 
least /i due to all intervening halos. 

The sum of all cross-sections in a fiat Universe can be 
written as 



dtot {n,Zs,Rs,ef,) ^ 47r 
dM 



^oj Jo 

a{fi, zg, Zs,M, Rs,e^)nc{zd, M)(l + z^fP'^iza) 

where nc{zd, M) is the number density of lenses of mass 
M at redshift Zd and D{zd) is the angular diameter dis- 
tance at redshift Zd- 
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Fig. 3. — Above curves demonstrate the correlation of finite 
source effects with halo ellipticity. As halos become more ellipti- 
cal they become more sensitive to finite source effects and dampen 
faster at high magnification. The colors correspond to different 
lens halo ellipticities: blue, cyan, green, and black correspond to 

= 0.1, = 0.2, = 0.3, and = 0.4 respectively. Solid lines 
are the cross-sections for a source with a radius of 0.023/?Bi„ while 
dashed lines correspond to a O.OAlREin radius source. 



Using the mass and redshift of our fiducial isothermal 
halo profile, to find the cross-section for each lens as a 
function of mass M and redshift Zd we can use the fol- 
lowing scaling relation 



a{^.,Zd,Zs,M, Rs.e^) 



(13) 
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where the subscript refers to the quantities for a nor- 
malization halo from which the cross-section for other 
halos are achieved by the above scaling relation. R'^ is 
the Einstein radius of the scaled halo and ctq is a look-up 
table of the cross-sections as a function of halo elliptic- 
ity and source-size to Einstein radius ratio. x(zd) is the 
scaling factor due to the dependence of the halo velocity 
dispersion, ct^ on redshift and is given by 



^ / HHzd) Ajzd) 
^^^"^ \m{zdfi) A(zrf,o) 



1/3 



(14) 



The probability that a source at Zs is magnified by a 
factor greater than /imm is then 



F(/i > ^ir,lin) = 



Cr{fJ, > flmin) 

47rL>? 



(15) 



Plotting dP/dzd for a range of expected source red- 
shifts is a good measure of the redshift distribution of 
lenses. Figure |4] shows that for sources beyond Zg — 1.0 
a significant population of lenses are located at redshifts 
of Zd ^ 0.4 to Zd ^ 1.0. 

5.4. Lensed Number Counts 
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Fig. 4. — Lensing probability dP/dz rescaled to 1 for visual com- 
parison of the peak. As sources move from Zs = 1.0 to much higher 
redshift of 2s = 5.0 the most dominant lens population rcdshift 
peak changes from ~ 0.4 to ~ 1.0. 

Gravitational lensing conserves the surface brightness 
of the lensed sources. Consequently any sources magni- 
fied by a factor of ^ are also times brighter. In addition 
to making sources appear brighter, gravitational lensing 
also dilutes the source populations by magnifying the ob- 
served solid angles. Therefore it also dilutes the number 
counts by a factor of ^. 

As discussed extensively in Jain and Lima (2009), these 
effects can be combined for a large survey to obtain the 
observed number counts as 



dn 
dS 



I dP dh ~ S 



■j,z)dzd^' (16) 



where we denoted the observed flux as S and the unlensed 
flux as S, such that S = fiS and the unlensed differential 
source sky density as dh/dS. 

The four panels in Figure [5] show the lensed count pre- 
dictions for SPT as the most relevant parameters are 
varied. The cosmological uncertainty is clearly subdom- 
inant to uncertainties in the source properties and lens 
properties. In particular, the overall normalization of 
the mass profiles (f^) and the source size can change the 
expected number of strong lenses by large amounts. 

Figure H] and [7] illustrate that there are lensing mod- 
els which fit the SPT source counts quite well for both 
the Durham and UBC unlensed source count models. 
The lensing parameters are slightly different, with the 
Durham model preferring slightly larger source sizes and 
smaller ellipticity. The contribution of various redshift 
ranges to the lensed number counts can also be seen to 
be slightly different for the two models. Figure [6] shows 
that the Durham model finds that most of the sources 
lie at z > 3 over the entire SPT flux range, while Figure 
[7] shows that the UBC model has a significant fraction of 
lensed sources as low as z ~ 2. 

In Figure [S] and [S] the distribution in magnification can 
be seen to be a strong function of flux for both models. 
At the low-flux end the cross-section is dominated by rel- 
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Fig. 5. — Effects of single parameter variations around a fiducial 
model (/o- = 1, (Tg = 0.801, Rs = ikpc, and e = 0.15) on predicted 
number counts of Durham unlensed model: cfi,[M) (top left), erg 
(top right), source size (bottom left) and ellipticity (bottom right). 

atively low magnifications (~ 10), while higher fluxes are 
increasingly dominated by larger magnifications. This is 
not surprising, given the steep unlensed luminosity func- 
tion. 

The mean magnification at different fiuxes is seen to 
have significant differences between the Durham and 
UBC models, suggesting that this could be a useful di- 
agnostic for reconstructing the unlensed luminosity func- 
tion. In Figure[5]the Durham model is largely dominated 
by the highest magnifications in the region of the SPT 
counts, while Figure |9] shows that the UBC model shows 
a transition from low magnification to high magnification 
in the flux range where SPT has reported constraints. 
This is a direct reflection of the differences in the un- 
lensed counts: the Durham model has an abrupt fall-off 
at the high flux end of the unlensed counts, well below 
the flux range probed by SPT, while the UBC model has 
more unlensed high redshift bright objects. 

6. DISCUSSION AND CONCLUSIONS 

The observed number counts of galaxies at mm wave- 
lengths are easily explained by gravitational lensing. 
However, the predictions for the number counts have sev- 
eral large sources of uncertainty, several of which we have 
explored in some detail. 

The physical size of sources is an important factor in 
determining the amount of magnification. This has been 
an im portant sources of u ncertainty in previous calcula- 
tions (jPaciga et al.l [20091 ) . Mathematical point sources 
can be treated analytically, but for any realistically ex- 
tended source the magnification can be strongly affected. 
For large sources (compared to the lens Einstein radius) 
the maximum magnification is considerably reduced. As 
the relevant scale is the Einstein radius (a property of the 
lens), this effect will depend strongly on the mass of the 
lens. At intermediate magnifications (i.e. 10 < /x < 40) 
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Fig. 6.— Best-fit lensed counts for SPT 220 GHz for Durham FiG. 8.— Magnification breakdown of the predicted SPT 220 

model. The lensing consists of halos with = 0.2 and Rs = 3kpc. GHz lensed counts for Durham model. 

The colored curves show the contribution of each redshift bin to 
the total counts. The unlensed counts only include sources with 
z > 0.2 and the data points are SPT 220 dust counts with IRAS 
counterparts removed. 
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Fig. 7. — Same as Fig. |6]for UBC model with best fit parameters 
= 0.3 and Rs = Ifcpc. 

finite source effects can increase tlic lensing cross-section 
because larger sources have higher probability of having 
a part of them fall within the inner caustics and be highly 
magnified. 

For the Durham model we observe that allowing 
sources with radius larger than 5 kpc leads to predicted 
counts falling below the observations. While it is not 
expected that the massively star-forming region in these 
sources is significantly larger than this, this is another ob- 
servational clue to the physical processes in these galax- 
ies. This has alrea dy been measured obs ervationally for 
individual sources |Momiian et al.l [20101 ) . More exten- 
sive future follow-up observation of more SMG lensed 
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Fig. 9.— Same as Fig. [8] for UBC model. 

sources could give better observational constraints on the 
range of SMG sizes. In addition it is possible to make a 
superposition of the curves of number counts for differ- 
ent assumed source sizes in order to achieve total counts 
for a diverse population of source sizes. This is mainly 
relevant to the distinction of AGN sources from star for- 
mation regions which extend over much larger spatial 
extents. For example assuming a 30% AGN contribution 
to the total counts one can superpose the curve of counts 
for point sources (AGN) and extended sources (star for- 
mation) with the weight of 0.3 and 0.7. Nevertheless 
this work shows that the effects of other uncertainties in 
modeling the lens (e.g., halo velocity dispersion normal- 
ization) on the counts dominate the more subtle effects 
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Fig. 10. — Normalized redshift distribution of the lensed sources 
at different observed flux levels. The model consists of Durham 
SPT 220 counts with z > 0.2 lensed with halos with an ellipticity 
of = 0.2 and sources with a radius of 3.0 kpc. 

such as the distribution of source sizes. 

EUipticity is important in two ways: it increases the 
total area of high magnification by increasing the area of 
the so-called "diamond caustic," but it also leads to more 
sharply-defined structures in the source plane, making 
finite source size effects more important. The trade-off 
between these effects was found to be an important effect, 
as this can change the local slope of the lensed number 
counts. 

A more realistic approach would be to assume a dis- 
tribution of ellipticities and source sizes to constrain the 
parameter space. However, with the strong parameter 
degeneracies we have identified it is clear that more data 
will be required to separate these effects. 

There remains significant uncertainty in the expected 
redshift distribution of the sources that comprise the 
mm-wave background from star-forming galaxies. Grav- 
itational lensing provides a magnified view of these ob- 
jects, providing an opportunity to probe the star forma- 
tion history. As a measure of the source redshift distri- 
bution we have plotted (Figure [T0|) the redshift distri- 
bution of the lensed counts (normalized to 1 for visual 
purposes) for three different observed flux levels. The 
redshift distribution of lensed sources at each observed 
flux level depends both on the redshift distribution of 
the unlensed model and the applied lensing probability 
function. 

Source redshift distribution is an important indicator 
of the nature of submm populati ons. The recent SPT ob- 
servations bv lVieira et al.l (|201Cl[ ) show a large population 
of mm sources which were previously unseen in other cat- 
alogs. These sources are speculated to be high-redshift 
gravitationally lensed galaxies. Our lensed number count 
plots suggest that this can be a possible scenario. Fig- 
ure [6] and [7] show that sources with redshifts larger than 
z > 3 (red curve) constitute at least half of the observe d 
non-IRAS-detected dusty sources of lVieira et all (IMoh . 

Throughout this work we have assumed SIS profiles 
as the mass model for lensing halos. Strong lensing is 



very sensitive to th e inner structure of the halos very 
close to the center (|Mead et al.ll2010( ). Rotation curves 
of galaxies which are approximately flat to large distances 
(| Rubin fc FordlHOTOl ) suggest that such a model is not a 
bad approximation to the mass density profile of galaxy- 
sized halos. While this is not likely to be a bad approxi- 
mation even up to cluster scales, the normalization of the 
mass density profile may not be expected to follow the 
self-similar scaling expected for the large-scale velocity 
dispersion of the dark halo. Given the strong sensitivity 
to the overall normalization of the density profile seen in 
Figure m strong lensing number counts will be sensitive 
to the details of the radial profiles of the mass density 
and its evolution with mass. 

A more re alistic calculation w ould involve taking an 
NFW halo (H avarro et al.l 119971 ). correcting for bary- 
onic condensation at the center, and populating it with 
small er subhalos a ccording to galaxy occupation num- 
bers (jOguril I2OO60 . The positioning of the subhalos 
could be an important factor, as would the density pro- 
files of the subst r uctur e. A similar work was done by 
iMeneghetti et all (|2003| ) where the effects of a cD galaxy 
on the cross-section of the main halo have been exam- 
ined. This could perhaps be better resolved by carrying 
out ray-tracing simulations through N-body simulations 
which include realistic baryonic mat t er. S uch a work 
has been carried out bv lHilbert et al.l ()2007l ) where they 
study ray-tracing through the Millennium simulation; 
they include baryonic matter based on semi-analytic 
mode ls of star formation in a related study (jHilbert et al.l 
120081 ). These studies have not had the required resolu- 
tion at galaxy-scale levels to study strong lensing due to 
substructure. In general, there is very little empirical 
guidance at this time for connecting Einstein radius to 
halo mass for halos smaller than galaxy clusters; popu- 
lation studies of strong lenses may provide some of the 
best constraints. 

Since the SPT does not have the power to resolve mul- 
tiple images of strongly lensed sources, the magnifica- 
tion of individual images does not affect the observa- 
tions. Hence in our calculations of magnification we have 
summed over the area of all the images of each source. 
If future observations succeed to resolve individual im- 
ages then the multiplicity of the images (e.g. the ratio 
of quadruply lensed to doubly lensed images) can put 
stronger constraints on the lens parameters, in particular 
lens ellipticity. The ratio of quadruply to doubly lensed 
sources can be obtained from the ratio o f the area of the 
diamond caustic to the radial caustic (jKormann et al.l 
Il994f ) which strongly depends on the ellipticity of the 
lens. There has however been a strong disagreement be- 
tween the observed ratio of t he quad to double image s 
and theoretical calculations (jRusin fc Tegmarklf200ll) . 
One of the possible effects causing this discrepancy may 
be observational selection effects due to the higher mag- 
nification of the quads. For example we calculate that 
the ratios of /i >= 10 cross-section (a typical magnifica- 
tion of SPT sources based on our work) to the area of 
the diamond caustic for ellipticities of e„ = 0.1, = 0.2 
and e«, = 0.3 are 4.12%, 19.24%, and 49.28% respectively. 
However perhaps the issue of quad to double ratios re- 
quires a more involved analysis which is beyond the scope 
of this work. 

The methods used here can be generally applied to 
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number count predictions for observations with other in- 
struments such as Herschel (operating at 250, 350, and 
500 ^m), BLAST, etc. The large sky coverage of SPT 
makes it particularly efficient, probing the rarest objects, 
where a large fraction are strongly lensed. On the other 
hand, the deeper Herschel observations allow counts to 
lower fluxes, allowing better characterization of the un- 
lensed source distribution. A combined analysis of Her- 
schel and SPT lensed counts will provide key insights into 
properties of both gravitational lenses and star-forming 
galaxies at high redshift. The results of iViei ra et al.l 
were based on less than 10% of the final pro- 
jected SPT survey area, and the Herschel observations 
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are just now starting to emerge. The upcoming Ata- 
cama Large Millimeter Array will have the resolution to 
make detailed images of these sources, providing more 
precise observational constraints on theoretical models 
presented here. 
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